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We present some basic elements of the treatment of particle multiplic- 
ities in jets from high energy collisions within perturbative QCD. Then 
we discuss the universal features of the inclusive particle spectrum for the 
limiting case of momentum p — > (or pt — > 0) as expected from soft QCD 
gluon bremsstrahlung. The energy independence of the invariant particle 
density in this limit Iq = E | p ^o is predicted as well as the dependence of 
this quantity on the the colour factors characteristic of the underlying par- 
tonic processes. These properties are hrst recalled from e + e~ collisions and 
then extended to pp and nuclear collisions according to Ref. pQ. Present 
data support these predictions. It will be interesting to see whether new 
incoherent contributions show up in the new energy regime of LHC. 

PACS numbers: 12.38 Bx, 12.38 Qk, 13.85 Hd 



1. Introduction 

The production of multi-hadron final states at high energies is described 
within QCD as a two-phase process: in the first phase there is some hard 
scattering (strong or electroweak) of the incoming elementary objects like 
quarks, leptons or gauge bosons. The produced quarks and gluons ("par- 
tons") will form jets of partons by gluon bremsstrahlung and quark pair 
production according to the rules of QCD perturbation theory for a char- 
acteristic cut-off scale Qq. In a second phase the partons reinteract and 
hadrons are formed which ultimately decay into stable particles. These pro- 
cesses are not accessible in perturbation theory and particular models are 
applied for their description. 

The simplest and best understood high energy process is e + e~ annihila- 
tion into hadrons. It is initiated by the process e + e~ —> qq which evolves into 
two hadronic jets dominantly. In pp collisions the protons in the primary 
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hard collision act as a collection of partons and in an event triggered for large 
transverse energy the partons scattered into large angles form sidewise jets 
while the spectator jets follow the direction of the incoming protons. The 
soft particles in this case form the "underlying event" and this phenomenon 
is under intense investigation today. In addition, there is the possibility 
of multiple independent parton-parton interactions considered important at 
the highest available energies. There are also the untriggered "minimum 
bias" events which may result from small angle parton-parton scatterings. 
Finally, in nucleus nucleus collisions there may occur hard parton parton 
interactions as in pp collisions, in addition there are parton reinter actions 
in the large nucleus, multiple nucleon-nucleon interactions and, with special 
interest, the new collective phenomena like quark gluon plasma formation. 

Although very different phenomena appear in the various processes there 
are some remarkable simplifications for very soft particles with momentum 
in the limit p — > in all processes. We consider for inclusive particle distri- 
butions the limit for the particle density 



E dN 
d?p 



(1) 



In this limit the Born term in the perturbative expansion dominates and 
this leads to some universal features for all processes 

1. inclusive spectra become energy independent 

2. the relative normalization of spectra in different processes is given by 
the colour factors relevant for the minimal partonic process. 

This holds for QCD partons, but we assume the same is true also for 
hadrons. 

These properties can be understood qualitatively as follows. A soft 
gluon is coherently emitted from all final partons. Having a large wave- 
length it cannot resolve any detailed intrinsic jet structure. It "sees" only 
the total colour charge which is carried by the primary partons, and these 
are represented by the Born term for the minimal partonic process in the 
perturbative expansion. 



2. Inclusive properties of QCD jets 

2.1. QCD evolution equations 

We begin by recalling the main tools to derive the inclusive observables 
for parton jets. They are obtained analytically in QCD using the concept 
of evolution equations (see, for example, Refs. [2J [3], some more recent 
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results will be added). Let us consider the partons in a jet emerging from 
a primary parton of energy E within the opening angle 0. First we con- 
sider the global observables like mean multiplicity (n), factorial moments 
fq = {n{n — 1) . . . (n — q + 1)) of the multiplicity distribution which can be 
derived from a generating function 

oo 

Z(Q,«)=£p n (Q)u n (2) 

n=l 

for the jet scale Q = EQ at small angles and the probability P n for 
production of n particles as 



dZ(Q,u) 

n = 



du 



_ d n Z(Q,u) 

Jq — 



(3) 

u=l 



In a corresponding way we obtain inclusive distributions D(k) = -pjjk i.e. 
the number of particles in the interval d 3 k, and, more generally, inclusive 
correlation functions D^ n \k\ . . . k n ) from a generating functional which de- 
pends on a probing function u(k) 

Z(Q, U ) = lZj d * kl ■ ■ ■ d3k n P n(h . . . k n )u{k X ) . . . U{k n ) (4) 



D^(ki ...k 



5 n Z({u}) 



nJ 5(u(k 1 )...5(u(k n ) 



(5) 

u=l 



where P n is the probability distribution of momenta /cj. 

For this generating function or functional an evolution equation is de- 
rived in the scale Q = EQ which yields, by appropriate differentiation, the 
equations for the observables like multiplicities and inclusive spectra. In 
differential form one finds the coupled equation for quark and gluon jets 
(a = <?, s) [U 2j0] They can be considered as extensions of the well known 
DGLAP evolution equations towards low particle energies taking into ac- 
count soft gluon coherence as realized in a probabilistic way by angular 
ordering [8] : 

1 — Zc 



— Z a (Y,u) = ^ / dz^—P bc (z) x 

6,c Zc 



{Z b (Y + In z, u)Z c {Y + ln(l - z),u) - Z a (Y, u)} (6) 
Z a (0,u) = u (7) 



1 Simplified forms have been obtained before. [B] 
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The evolution variable is taken as Y = ln(E@ /Qq) with the non-perturbative 
hf cut-off Qq; the argument of the running coupling is kx = min(z, 1— z)EQ. 
The evolution equation (J6j) describes the decay of a parton jet a at scale E 
into two parton jets b, c at scale zE and (1— z)E with probability Pbc(z), the 
so-called DGLAP splitting functions. The second equation ([7]) represents 
the initial conditions at threshold (EQ = Qq) and means that the parton 
jets start just with the initial parton a. 

Asymptotic solutions can be obtained in the Double Logarithmic Ap- 
proximation (DLA) which includes only the dominant contributions from 
the singularities at small angles and energies in the emission probability. In 
this approximation the splitting function P ga (z) ~ 1/z in (|6|); the next to 
leading single logarithmic terms are included in the Modified Leading Log- 
arithmic Approximation (MLLA). Up to this order the results from Eq. ([B]) 
are complete; further logarithmic contributions beyond MLLA can be calcu- 
lated, but they are not complete and neglect in particular process dependent 
large angle emissions. Nevertheless they improve the results considerably as 
they take into account energy conservation with increasing accuracy. The 
full solution of Eq. ([6]), corresponding to the summation of all logarithmic 
orders can be obtained numerically. Alternatively, one may calculate results 
of the QCD cascade from a Monte Carlo generator, such as ARIADNE [9j, 
which applies the same kx cut-off procedure as Eq. Q. 

2.2. Parton Hadron Duality Approaches 

So far we have discussed the properties of a jet of partons obtained from 
perturbation theory using an artificial cut-off at low scales Qq. The appli- 
cation to multiparticle observables needs an additional assumption about 
the hadronization process at large distances which is governed by the color- 
confinement forces not accessible by perturbation theory. 

The simplest idea is to treat hadronization as long-distance process, 
involving only small momentum transfers, and to compare directly the per- 
turbative predictions at the partonic level with the corresponding measure- 
ments at the hadronic level. This can be applied at first to the total cross 
sections, where at the low energies the resonance structures are represented 
in an average sense. The perturbative approach also describes jet produc- 
tion for a given resolution; here the collection of partons is compared to 
hadronic jets at the same resolution and kinematics. This approach has 
led to spectacular successes and has built up our present confidence in the 
correctness of QCD as the theory of strong interactions. 

In a next step one may carry on such a dual correspondence further to 
the level of partons and final hadrons. This procedure turns out successful 
for "infrared and collinear safe" observables which do not change if a soft 
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particle is added or one particle splits into two collinear particles. Such 
observables become insensitive to the cut-off Qo for small Qq. Quantities of 
this type are energy-flows and -correlations as well as global event shapes 
like thrust etc. 

Further on, one may compare partons and hadrons for observables which 
count individual particles, for example, particle multiplicities, inclusive spec- 
tra and multiparton correlations. Such observables depend explicitly on the 
cut-off Qo (the smaller the cut-off, the larger the particle multiplicity). 

According to the hypothesis of Local Parton Hadron Duality (LPHD) [TU] 
the hadron spectra are proportional to the parton spectra where the con- 
version of partons into hadrons occurs at a low virtuality scale, of the order 
of hadronic masses, i.e. Qo ~ f ew hundred MeV, independent of the scale 
of the primary hard process. 

While this hypothesis has been suggested originally for single inclusive 
spectra it can be generalized to more complex situations of the form [3] 

0{x x , X2 , ■ ■ ■ ) | hadrons 

K 0(xi,X 2 , - ■ ■ ,Qo,A)\partons, (8) 

where the non-perturbative cut-off Qo and the "conversion coefficient" K 
have to be determined by experiment. The conversion coefficient should be 
a true constant independent of the hardness of the underlying process. 

In a more recent analysis mean multiplicities and higher multiplicity 
moments have been calculated both for sub-jets of variable cut-off scale 
Q c ("jet virtuality") and for hadrons with cut-off Qo m e + e~ annihilation 
[1H \I2\ with the smooth transition from jets to hadrons for Q c — > Qo0 
For jets at fixed cut-off Q c the normalisation is K = 1 in ((5]). It turns out 
that a unified description of jets and hadrons was possible with the common 
normalization 

Km I. (9) 

In this case the hadronic cascade has been represented by the partonic 
cascade in the average with the same multiplicity of partons and hadrons. 
So the resonance bumps in the multi-particle spectra are just represented by 
the corresponding smooth perturbative spectra in the average. This parton- 
hadron-jet correspondence implies that a hadron corresponds to a parton jet 
of resolution Qq. 

When comparing differential parton and hadron distributions there can 
be a mismatch near the soft limit caused by the mass effects (partons are 
taken as massless in general). This mismatch can be avoided by a proper 
choice of energy and momentum variables. In a simple model [THE] partons 
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and hadrons are compared at the same energy (or transverse mass) using 
an effective mass Qo for the hadrons, i.e. 



ET,parton — hr,parton ^ Ex^hadron — \J ^T,hadron + Qo> U^) 

then, the corresponding lower limits are kT, pa rton Qo and kx^hadron 0. 

We should remark that these duality approaches are justified primarly 
by their phenomenological success and their intrinsic simplicity and not yet 
by a convincing theoretical derivation from QCD. In particular they allow 
compact analytical solutions for the observable quantities in the available 
approximations to QCD (DLA, MLLA. ..), which is not possible for the 
phenomenological hadronization models of high complexity. 



3. Quark and gluon jets: global observables vs. soft limit 

3.1. Global observable: particle multiplicities 

Before we derive the soft properties of particle spectra we discuss the 
mean particle multiplicity in a jet as the most simple example of a global 
event characteristic. Here the higher orders in the QCD perturbation theory 
are very important. By differentiation (|3j) of the evolution equations of the 
generating functions (|6j) one obtaines the evolution equations for the parton 
multiplicities J\f a in quark and gluon jets (a = q, g) [TT| 

dM a {Y) 1 ^ f 1 , a s (k T ) 



dY 

b,c 

x [M b (Y + \nz)+N c {Y + \n{l-z))-N a {Y)} (11) 
with initial conditions 

N a (Y)\ Y=Q = 1, (12) 

which imply there is only one particle in a jet at threshold. Starting with 
this initial condition one can obtain by iteration of the evolution equation 
the perturbative expansion. 

The asymptotic behaviour can be derived from (jlip using the ansatz 



A/- 9 (Y)~exp( / \(y)dy), (13) 



where the anomalous dimension 7 has an expansion in 70 = ^J2Ncoi s /'k 

7 = 7o(l - ai7o - a 2 7o - °37o • • •)> ( 14 ) 
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C A /C F = 9/4 
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Fig. 1. The ratio of the mean multiplicities in gluon jets and quark jets N g and N q 
obtained from e + e~ experiments; results from perturbative QCD show the large 
higher order corrections for a global observable (from [16]) . 



and for the ratio of gluon and quark jet multiplicities 

_Mg _C A 



Nq C F 



(1 - ri7o - r 2 7 - r 3 -y . . .) 



(15) 



with QCD colour factors 



C A = N C = 3, C F = -. 



(16) 



The coefficients dj and n can be derived from the evolution equations. At 
high energies the leading behaviour in MLLA for both quark and gluon jet 
multiplicities is given by 



\nM(Y) ~ c 1 / V / a s (Y) + c 2 In a s (Y) + c 3 
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C-2 



1 10 /L 

4 + 27^ /6 ' 



C A 



-n 



(17) 
(18) 



with the arbitrary constant C3, and this behaviour describes well the data 
in e + e~ annihilation at LEP-1 and LEP-2, for review, see |17j . more recent 
results have been presented by DELPHI [ig [H [20] and OPAL [2H[22]. 

The important role of higher logarithmic orders can be studied in the 
behaviour of the multiplicity ratio r in (|15|) . The asymptotic limit r = 



Ca/Cf acquires large finite energy corrections in NLLO [23, 24J and 2NLLO 
order [251 [26] 



with hi = ^, h 2 = §2 r, h 3 = and ^4 = T' also 3NLLO results 

have been derived [27]. Results from these approximations [T7J are shown 
in Fig. 1 together with the numerical solution of the MLLA evolution equa- 
tions (jlip obtained in 1998 [llj . which takes into account all higher order 
corrections from this equation and fulfils the (non-perturbative) boundary 
condition (|12p . All curves are absolute predictions, as the parameter A (and 
Qo in case of the numerical calculation) is adjusted from the growth of the 
total particle multiplicity in the e + e~ jets. The slow convergence of this 
y / aj expansion can be seen and there are still considerable effects beyond 
3NLLO. The numerical solution is also in close agreement with the MC 
result at the parton level obtained [22J from the HERWIG MC above the 
jet energy £j et > 15 GeV (£j e t = Q/2 in e + e~ annihilation) and ~ 20% 
larger at E- ]et ~ 5 GeV. This overall agreement suggests that the effects not 
included in the MLLA evolution equation, such as large angle emission, are 
small. 

These numerical results are also compared in Fig. 1 with data from 
OPAL [22] where the data on gluon jets are derived from 3-jet events in 
e + e _ -annihilation. Note also that a proportionality constant K relating 
partons and hadrons according to LPHD drops in the ratio r. The results 
obtained from DELPHI [TjJ] fall slightly below the curve by about 20% at 
the lowest energies but converge for the higher ones; the CDF collaboration 
comparing quark and gluon jets at high in pp collisions [28J finds the 
ratio r in the range 5 < Ejet < 15 GeV a bit larger, closer to the 3NLLO 
prediction, but with larger errors and therefore still consistent with the LEP 
results. 



Next, we consider the inclusive distribution D(£,Y) of partons in the 
momentum fraction x = k/E or £ = ln(l/x) within a jet with primary 
parton energy E and opening angle 0. The evolution equation for D can 
be obtained by functional differentiation of ([6]) (for a review, see |3jb). At 
small x (large £) the angular ordering [7j of the cascade evolution which 
takes into account the soft gluon interferences in a probabilistic way plays 




(19) 



(20) 



3.2. Inclusive energy spectrum: soft limit 
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an important role. For large x the equations approach the well known 
DGLAP evolution equations. 

For the present discussion we restrict ourselves to the simplest approxi- 
mation, the DLA with fixed coupling, where only the most singular terms in 
the splitting functions for g — > gg and q —> qg are kept. Then the evolution 



equation for parton a reads 

D^,Y) = 8*M) + f V [ C dy^ 7 2 (y)Dl(Z',y). (21) 



This equation can be solved by iteration. For fixed 70 ~ y/cts one obtains 
the perturbative expansion 

Bl&Y) = mi) + ^ll{Y-0 + \^lfc(Y-i? + --- (22) 

Ca \Y-Z 



= mo + ^7 y /i(2 7 oVe(^-o) (23) 

with the modified Bessel function I\. The coherent soft gluon emission 
leads to a depletion of the spectrum at large £, also called "the hump- 
backed plateau". From the inclusive distribution (|22p one can obtain the 
double differential distribution in energy and angle by the differentiation 
over Y which yields = + %^lt^(Y ~ + • • • > or ! in tne original 

variables 

dN a 2C a AN c C a 2 E k T 

^@ = as + — k@ a ° ln I ln q~ + • ■ ■ (24) 

Here we recognize in the leading term of 0(a s ) the well known Born term 
for soft gluon bremsstrahlung as in QED, but with the appropriate QCD 
colour factors. One observes that only the Born term survives in the soft 
limit where kx — > Qo; in this limit we find the simple universal properties 
emphasised in the introduction: the particle density becomes independent 
of energy E and is proportional to the relevant colour factor for the minimal 
process, that is here the gluon emission from the quark or gluon jet with 
C a = Cf or C a = Ca respectively. 

This result can be generalised to the accuracy of DLA with running 
coupling which is obtained by iterating ()2ip accordingly up to 0(a 2 s ) which 
is appropriate for the low momentum region, furthermore results within 
MLLA have been derived as well |15| . The above properties in the soft limit 
remain unaltered. A comparison of these calculations is shown in Fig. [2] 
where a model dependent kinematic relation as in (jlOp is used. One can 
see that the data are rather well described by the model which predicts an 
energy independent particle density in the soft limit p — > 0. 
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Fig. 2. Inclusive particle densities 



' :i /< 



v^[GeV] 

of charged particles in e + e~ annihila- 



tions; left panel: momentum spectra for different energies, right panel energy ^fs- 
dependence for fixed momenta p, lower right panel: kaons and protons; calculations 
within MLLA accuracy. 



3. 3. Colour factors in quark and gluon jets 

In order to check the sensitivity to the colour factors in the Born term we 
should study the dependence of the soft limit of momentum spectra in quark 
and gluon jets. Separating quark and gluon jets results in uncertainties 
increasing with smaller momenta. Therefore an alternative procedure has 
been suggested [15] which studies the radiation into a cone perpendicular 
to the production plane of 3-jet events in e + e~annihilation. 

One can consider two extreme limits with two jets aligned: 

a) a quark and a gluon jet are parallel and recoil against a quark jet, in 
this configuration the soft perpendicular radiation cannot separate the two 
parallel jets and the intensity is as in a qq dipole, proportional to the colour 
factor Cf; 

b) the quark and antiquark are parallel and recoil against the gluon, in this 
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case the primary configuration acts like a gg dipole with colour factor Ca- 
Of course, in a realistic experiment one cannot go to such extreme limits 
because of the finite width of the jets, but one can measure a certain range 
of the angles in between the jets which interpolates between the limits. 

The soft gluon bremsstrahlung from a "colour dipole antenna" {qq or 
gg) is given by [2] 



=(£r% W ^. W'..H" S ) = 2C,. F (Q) (25) 

(Q) = J, ''Z 6 " n (26) 

(1 — cos e>t S )(l — cos Wj s J 

where 0y is the angle between partons i and j and s denotes the soft gluon. 
Then, for the aligned qq antenna one obtains Wf = ACf/ sin 2 Q qs . 
In a 3-jet event e + e~ —> qqg one finds in lowest order 

w qqg {n g ) = c A l(€g) + ifg) - ^ ( 27 ) 

i.e. there are two dipoles between each of the quarks and the gluon of 
strength Ca and a colour suppressed dipole between the q and the q. For 
q\\g one finds W = ACpj sin 2 Q qs like a qq dipole, and for q\\q one obtains 
W = AC 'a/ sin 2 Q gs like a gg dipole, as anticipated above. 

The radiation perpendicular to the 3-jet plane (cos 0j s = 0) normalised 
to the same radiation in 2-jet events is then given by the simple formula 

^£ - ET = ^ r(e ..) (28) 
r(@ij) = ^[(l-cose g5 ) + (l-cose 59 )-^(l-cose gg -)] (29) 

Such a measurement has been carried out by the DELPHI collaboration 
[29} 130] . The above formulae should apply in the soft limit where the Born 
term of 0(a s ) dominates. The pt spectra in the cone perpendicular to 
the 3-jet plane are found all very similar for px ^ 1 GeV, therefore one 
can study instead the integrated multiplicity in the respective cones. One 
observes first that the multiplicities in the cone is well described by the 
above formula (|29p . the data are accurate enough to even notice the 1/Nq 
term in (|29p . Furthermore, the data for the ratio on the l.h.s. of (|28p are 
found to depend linearly on the function r{®ij) and one obtains from the 
slope 

= 2.211 ± 0.0U(stat.) ± 0m3(syst.) (30) 

Up 



12 



which is well consistent with the expected Ca/Cf = 9/4 in QCD. 

From these studies of jets in e + e~annihilation we can conclude that the 
soft particle density indeed follows the prediction of the soft gluon Born 
terms emphasized in the introduction 

a) The spectra become independent of energy for p — > 0; 

b) the soft particle density varies with the orientation of the colour an- 
tenna as predicted, this implies that the soft particle density in quark and 
gluon jets becomes proportional to the colour factors Ca and Cp. This 
is in strong contrast to the behaviour of global observables like the mean 
multiplicity, which obtains large higher order corrections from the integral 
over the perturbative expansion such as (|24p and the ratio r = N g /N q is as 
low as 1.5 at LEP energies instead of 2.25 (see Fig. 1). In the soft limit 
there is no phase space for subsequent emissions, nor for energy momentum 
conservation effects. With this experience we now investigate the hadronic 
collisions. 

4. High energy pp collisions 

We discuss here the "minimum bias" events, which we consider as non- 
diffractive events. In order to estimate the very soft particle production we 
look for the minimal partonic process which could be responsible for the soft 
gluon bremsstrahlung. We assume that the relevant process is a semihard 
2 — > 2 + g s scattering of lowest perturbative order where any two partons 
inside the proton can scatter with one- gluon exchange at small angles. The 
exchange of a colour octet gluon at small angle transfers the colour from the 
colour singlet protons to the two outgoing partonic clusters which are the 
colour octet sources of soft gluon bremsstrahlung. In the minimal configu- 
ration each proton splits into a quark-diquark pair which scatter by gluon 
exchange. In large Nq approximation the process corresponds to two radi- 
ating colour triplet antennae responsible for bremsstrahlung from initial and 
final partons. It should be added that also more complex partonic processes 
will end up in the production of two colour octet systems as discussed in [1] . 

Therefore, according to our general rules, we expect the energy inde- 
pendent limiting soft radiation density in pp collisions Iq P for p — > and, 
furthermore, we expect this density to occur in a fixed ratio to the corre- 
sponding density in e + e~ collisions as ratio of colour octet and colour triplet 
dipole sources 

p^O: I P P /Io +e '^ C A /C F , (31) 

just like the ratio of the spectra in gluon and quark jets discussed in the 
last section. 

This kind of relation (|3ip has been suggested by Brodsky and Gunion 
already in 1976 [3"T], but relating the integrated multiplicities in the central 
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rapidity region to these colour factors. From our QCD analysis we find 
Eq. (|3ip to be valid only in the soft limit while the ratio of integrated 
multiplicities is found closer to unity (see below). Similarly, relations of the 
kind (|3ip appear in some early phenomenological models, but again for the 
integrated densities only, as outlined in pQ. 

Our expectation of an energy independent Iq P is based on a coherent 
process. It would be violated if there were multiple parton-parton interac- 
tions (processes like 4 — > 4 + g s ) added incoherently. Such processes appear 
in some models at high energies (see, for example, PYTHIA [32]) and so 
the measurements at LHC can shed some light onto the contributions from 
such processes. 

We have studied the energy dependence of I pp using the results of 
fits to the invariant cross sections E-jjj^ measured in the energy range 
y/s = 20 . . . 1800 GeV obtained from the colliders at CERN, Fermilab and 
Brookhaven and measurements of inelastic cross sections Ui n . The pr spec- 
tra look qualitatively similar to those in Fig. [2] for e + e~ annihilations con- 
verging towards small pt but falling more steeply at high momentum. The 
Pt spectra have been fitted to distributions which at small pt behave like 



These fits are good down to the smallest measured pr ~ 0.1 GeV. Then 
one finds I pp = A/<Ji n from extrapolation pt —> 0. The functional form 
(|32|) is not analytic at px = and is therefore theoretically not satisfactory. 
This problem is avoided using the "thermal" parametrisation in terms of 
jut instead of pr 



as applied by PHOBOS [33] in nuclear collisions and a good fit down to the 
smaller pt ~ 0.03 GeV has been obtained. The extrapolated values Iq P are 
smaller by about 25% as compared to the fit ([32]) . 

The results from the available published exponential extrapolations are 
shown in Fig. [3l One observes a rather flat energy dependence over the two 
decades in energy with an average I pp ~ (7 db 1) GeV~ 2 . Note, that over 
this energy range the rapidity density ^~ would rise by about a factor 2. 
The different extrapolations f|32[) or f|33|) should not affect the trend of the 
energy dependence. For comparison with e + e~ annihilation it is better to 
normalise by the non-diffractive cross section which we take as 15% lower 
than the inelastic one, which yields ^Qnd ~ ± 1) GeV -2 . Then for the 



E—r- = Aexp(Bp T + ...)■ 



(32) 




(33) 
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Fig. 3. Soft limit of the invariant density E-^ of charged particles [(h + + h )/2] 
in pp collisions as a function of c.m.s. energy y/s (from exponential extrapolation). 



thermal fit and non-diffractive (minimum bias) events we obtain 

if« (6±1) GeV" 2 . (34) 

So far, there is not yet a fit result from LHC to be used for comparison. 

Next, we also compare this result with the soft limit in e + e~ annihilation 
to test our prediction (|3T]h We present two results. 

1. There is a dedicated comparison by the TPC collaboration |34j who 
compared their own data on e + e~ annihilation with pp data from the British 
Scandinavian collaboration on the pt spectra of the invariant density. 
The TPC data are presented as function of pr as determined from the 
sphericity jet axis. In the average over px both data sets for pions, kaons 
and protons are similar. A closer look, however, reveals, that the spectra 
fall more steeply with px in the pp collisions and there is a cross over of 
the spectra at low px- The appropriate extrapolation down to px near zero 
yields a larger density for pp collisions by a factor 2.0 — 2.7 depending on 
the kind of fit. 

2. Most other experiments present fits to the spectra in particle energy 
E (not pt)- Using the fit results from various experiments in the range 
yfs = 10 ... 29 GeV yields lf e ~ » (3.3 ± 0.5) GeV" 2 or the ratio 

I^/lf e ~ « (1.8 ± 0.4) -=- (2.4 ± 0.5), (35) 
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where the first (preferred) number refers to the thermal and the second to 
the exponential extrapolation. This result is consistent with our QCD based 
expectation for this ratio Ca/Cf = 2.25. 



5. Nucleus-nucleus scattering 

For the nucleus-nucleus (AA) cross sections we may consider two limiting 
cases in the relation to the pp cross section. 

1. In case of a point like interaction the particle densities in nuclear 
collisions are obtained by rescaling the densities in pp collisions by N co n, the 
number of nucleon-nucleon collisions, or, "the nuclear modification factor" 

pJVcoii = 1 dN AA /dp T . . 

Uaa N coll dN pp /dp T ' W 

is unity. The number N co u can be obtained from the Glauber model. 

2. In case of soft particle production we expect that such particles with 
a large wave length 1/pr ^ r are coherently emitted over a range r (from 
a nucleon rjy ~ l/ m 7r or a nucleus r A ~ 1/(30 MeV), which results in 
a reduced rate. Indeed, the inspection of the RHIC data [33j 133 EH [3H] 
shows the ratio R AA " falling below unity for small pr- An alternative way 
presenting data in the soft region is the normalisation to the number of 
"participating nucleons" 



_1 dN AA /dp 7 

(N part /2) dN pp /dp T 



= /AT 1 /n , (37) 



This concept has been introduced already in 1976 by Bialas, Bleszynski 
and Czyz [39], who found the number of participating nucleons, called 
there "wounded nucleons" as relevant scaling factor for soft production 
(fi^° rI « 1), i.e. each interacting nucleon should be counted only once 
and the rescatterings of the same nucleon be disregarded. 

Again we consider the energy dependence and the normalization of the 
particle production at central rapidity in the limit pt — > 0. A detailed 
study for various centralities (from peripheral to central AA collisions) by 
PHOBOS [33] shows that the inclusive pt spectra in normalisation (|36p or 
(|37p approach about the same densities at 200 and at 62.4 GeV. This implies 
that the energy dependence is the same for nuclei and protons, i.e. there is 
no sizable energy dependence. 

Concerning the normalisation, the same data [33] show that the quantity 
(|37p approaches unity for all centralities (within about 30%) at the lower 
limit of pt ~ 200 MeV. The STAR collaboration [3D] has measured this 
quantity with high precision down to 500 MeV and an extrapolation to the 
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"participant scaling" for px = is indicated. We combined fits of PHOBOS 
low pt and STAR AuAu data at 200 GeV using the thermal parametrisation 
and also the STAR pp data which yields the result 

l£ A /ll P ~ 160 ± 17, (38) 

which agrees with the calculated N part /2 = 172 (±15%), and is therefore 
consistent with 

p T ^0: <^^1 and l£ uAu « ^ I™. (39) 

This density for AuAu collisions is about six times smaller than expected 
for an incoherent superposition of collisions with N co u = 1040, where this 
number is obtained from Glauber model calculations. It is remarkable, 
that the "wounded nucleon" model works to the precision of about 10%, 
the accuracy of measurements and theoretical calculations. Note that this 
agreement is obtained only in the limit pt — > as can be seen from the 
STAR data [ID] : already for px = 0.5 GeV the deviation from "participant 
scaling" amounts to about 50%. 

How can this scaling result be understood? An incoming nucleon scat- 
ters successively at a number of nucleons in the nucleus (see Fig. U] b,c for 
representative diagrams of pA scattering) . The successive gluon exchanges 
yield again an outgoing colour octet state as in pp scattering, such that the 
rescatterings of the nucleon are not causing any production of additional 
particles. This happens if in the low px interaction only a quark and the di- 
quark appear as active partons, so that also a multi-gluon exchange cannot 
produce a higher colour multiplet than an octet H Alternatively, one can 
think of a larger number of exchanged gluons but taking into account that 
the colour octet exchange gives the dominant Leading Logarithmic contribu- 
tion both in DGLAP and BFKL kinematics. Then, from both viewpoints, 
each scattered nucleon produces dominantly a colour octet flow as in pp 
interactions, in agreement with the phenomenological result (|39|) . 

6. Universal composition of softly produced particles 

Finally, we may ask, whether the universal production of the soft parti- 
cles from gluon bremsstrahlung also reflects in their composition as detected 
by particle ratios. Such a universality can be expected if the source of the 
bremsstrahlung are the colour triplet dipoles generated in pp and AA colli- 
sions by gluon exchange. In that case not only hadronic collisions but also 
e + e _ annihilations have universal dipole sources. 

3 A model based on wounded quarks and diquarks has been developed in Ref. [3T], but 
for the description of pr-integrated rapidity distributions. 
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(b) 




Fig. 4. Diagrams, contributing to the pp and pA collisions in the minimal model 
for soft particle production: (a) In the pp collisions the exchanged gluon interacts 
with the colour triplet constituents q or qq in the proton to form an outgoing 
colour octet system; (b,c) In the pA collisions the proton can rescatter inside the 
nucleus and then forms a colour octet system again. This implies that the multiple 
gluon exchange acts as a single gluon exchange in the particle production. In this 
example: N coU = 3, N part /2 = 2. 

A similarity of particle ratios K/tt and p/ix in the e + e~ and pp reactions 
at pt < 0.5 GeV has been indeed noted already some time ago by the TPC 
collaboration [34 . In this measurement the transverse momentum pj- for 
the particle collisions was defined with respect to the sphericity axis. 

The pt dependence of the particle ratios for several hadronic collisions 
have been compared by PHENIX [32]. While at the large pr > 2 GeV the 
ratios p/n and K/tt tend to approach large values ~ 1 in the central AuAu 
collisions, these ratios are reduced for non-central and minimum bias pp 
collisions. Remarkably, these ratios converge for all the different processes 
towards lower px < 1 GeV. In Fig. [5] we collect data in the low px region on 
the ratio K~ /ir~ from the e + e~, pp and AA interactions. As one can see, 
these particle ratios, indeed, approach each other towards low pt < 0.4 GeV 
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Fig. 5. Convergence of particle ratios K jn towards small px for various pro- 
cesses: e + e~ annihilation (TPC [3D data, pr with respect to sphericity axis), pp 
(minimum bias) and central (0-5%) AuAu collisions (PHENIX 42 and STAR [43] 
Collaborations). 



pointing towards a dominance of multiple qq dipole radiation in all processes. 



7. Summary 

We note some universal features of the particle production in the limit 
P,Pt which we derive from the dominance of coherent QCD gluon 
bremsstrahlung in this limit. We consider the particle density Jo in this 
limit for which we predict 

1) the energy independence, and 

2) the dependence on colour factors according to the minimal partonic 
process (Born-term): 



a) e + e annihilation: Iq r 



lll^ = C A /C F ; 



b) pp scattering: Iq P /Iq e = Ca/Cf] 

c) AA scattering: I AA /I% P = (N part /2) C A /C A . 

These expectations are well met by the data. In consequence, the soft 
particles do not follow a universal thermal behaviour independent of the 
initial state. 

There is also some universality in the particle ratios which tend to con- 
verge to those from qq dipoles. Soft hadrons in the central region are pro- 
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duced first. In A A collisions these slow particles stay behind and do not 
participate in the equilibration. 

It will be interesting to study the soft limit at LHC energies. If there 
are new incoherent sources, as expected in some models with multiple in- 
teractions, the soft density Iq could start rising with energy. 
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